A Nano-Quantum Photonic Model for Justification of Dispersion in 

Single Crystal Film of NPP 



o 
o 

(N 
> 

o 

00 
(N 

(N 
> 

in 

(N 



o 



Hassan Kaatuzian, AliAkbar Wahedy Zarch * 
Photonics Research Laboratory 
Electrical Engineering Department 
AmirKabir University of Technology, Tehran, IRAN 
hsnkat o@ aut . ac . ir 
fat emahali@noavar . com 

February 1, 2008 



Abstract 

In this paper, we present a nano-quantum photonic 
model for justification of normal dispersion in a thin 
crystal film of NPP. In this method, we assume a laser 
beam consists of a flow of energetic particles. By precise 
analyzing of photon interaction with 7r-electron system 
of benzene ring in NPP crystal, we will attain refractive 
index (RI) in any wavelength and compare the results 
with experimental data. 
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1 Introduction 

Organic optical materials like MNA, NPP, MAP have a 
high figure of merit in optical properties in comparison 
with inorganic optical materials such as LiNbOs, ADP, 
KDP and GaP [Ij. Additionally N-(4-nitrophenyl)-L- 
prolinol (NPP) after MNA has the highest figure of merit 
between organic nonlinear optical materials [2j. Thus it 
is used for electro-optic and nonlinear optic applications 
[IlEllIlSllllElHElinillOlIIB. It may say that refrac- 
tive index (RI) is the most effective element in optical 
phenomena. Ledoux et al. first proposed a Sellmeier set 
for RI of NPP, based on measured refractive indices [3j. 
Banfi, Datta and co-workers design a more accurate Sell- 
meier set for RI of NPP for nonlinear optical studies 
[H |5l 6J; these data based on classic and macroscopic 
measurements and approaches. Some authors explain 
RI in molecular bases [TH |l3l |14] , although no data on 
real material exist in this literature. Our purpose in this 
paper is the explanation of RI or dispersion effect of a 
real material, a single crystal film of NPP, in an approach 
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that name nano-quantum photonic [TF^ IT6l [17] . This ap- 
proach is based on four elements: 1- quasi-quantum prin- 
ciple for justification of phenomenon in molecular scales, 

2- Knowledge of crystal network and its space shape, 

3- Short range intramolecular and intermolecular forces, 

4- Monte-Carlo time domain simulation. We suppose a 
laser beam is a fiow of photons when passes through sin- 
gle crystal film, interacts with delocalization 7r-electron 
system of NPP molecule and delays the photon in ev- 
ery layer. By precise calculation of these retardation in 
every layer, we obtain refractive index (RI) in specific 
wavelength. The results obtained from this method are 
agreeable to experimental data. Our previous work was 
about MNA aromatic crystal, with a simpler structure 
than NPP molecule, [iTl [18]. By simulation on more 
complex crystal, NPP, and considering the agreement 
of simulation results with experimental and Sellmeier 
data, the validity of assumptions and simulations will be 
proved. Hence, first we introduce crystal and molecular 
structure of NPP and comparison to Benzene molecule, 
after which we explain normal dispersion formula ac- 
cording to classical physics, electromagnetism, quantum 
electrodynamics and our microscopic model. In fourth 
section, we simulate and calculate the probability den- 
sity function (PDF) of electron presence in benzene ring 
of NPP molecule, approximately. In fifth section, we ex- 
press detailed calculation of RI. In sixth part, we clarify 
how to obtain RI for NPP crystal and compare simula- 
tion results with experimental and Sellmeier data. 



2 Crystal and Molecular struc- 
ture of NPP 

Organic molecular units and conjugated polymer chains 
possessing 7r-electron systems usually form as cen- 
trosymmetric structures and thus, in the electric dipole 
approximation, would not show any nonlinear second 
order optical properties. The necessary acentric may be 
provided by first distorting the 7r-electron system by in- 




Figure 3: Electron cloud for Benzene molecule that obtained 




Figure 1: The molecular compound of NPP. 



teraction with strong electron donor and acceptor groups 
[19] . In NPP molecule nitro group acts as an acceptor 
and prolinol group on the other side of benzene ring acts 
as a week donor (see fig. 1). NPP (Ciii^i4iV203) (Fig.l) 
crystallizes in the solid state in an acentric monoclinic 
(with space group P2i) structure and their parame- 
ters are:a=5.261A°,b=14.908A°,c=7.185A° ,/3=105.18°. 
The volume of unit cell is V=543.8(74° ). The molec- 
ular weight of NPP is M=222 and the crystal with 
two molecules in the unit cell has a calculated density 
D=1.36gm.cm-^ [Ij. The melting point of NPP is 116°^ 
and in the wavelength range of 0.48 to 2/i is transpar- 
ent pj. The most interesting property of NPP crystal 
is the proximity of the mean plane of molecule with the 
crystallographic plane (101); the angle between both of 
these planes being 11°. Nitro group of one molecule in 
downward connects to prolinol group in upper by hydro- 
gen bonding. The angle between b orientation of crystal 
and N(l)-N(2) axis is equal to 58.6° [Ij. Fig. 2 shows 
the crystal packing of the NPP. For accurate and valid 
simulation, these properties and angles have to be ex- 
erted. For benzene molecule, benzene ring is a circle 
(see fig. 3); Fig. 4 demonstrates the bond lengths of 
the NPP molecule. As we see in this figure the bond 
lengths in benzene ring are not same. In our simulation, 
for similarity we consider an ellipse correspond to circle 
for electron cloud. We obtained £=0.26 for that ellipse 
from simulation, fig. 5 shows this comparison. 



OO 

Huckle theory [20j and its approximation. 




Figure 2: The crystal packing of NPP ^. 




Figure 4: The bond lengths and angles of NPP molecule 
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Quantum Electrodynamics attains a famous disper- 
sion relation in form of: 



heo 
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by Density Matrix Method where: 
where W is the energy level and 



Mifc — Pik 
that introduce the notation: 



0.2 0.4 0.6 



for unit of the polarization of the incident and scattered 

Figure 5: The comparison between circles for Ben- 



zene molecule electron cloud and NPP molecule electron 
cloud, (approximately). 

3 Classic and microscopic model 
for dispersion 

Cauchy suggested an equation that performs a relation 
between wavelength and RI in form [21j : 



1 = 



A B 
A2 ^ 



(1) 



If we have three RI in three wavelength according to 
experimental data, then we have RI in any wavelength 
in transparent range of that substance. Ledoux et.al 
[3j proposed a Sellmeier formula for NPP crystal in the 
form: 
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and Datta et al. 
in the form: 



5 demonstrated a more accurate of it 



n' = A- 



B 



D 



1 



c 



1 



E 



that A,B,C,D and E is exerted from table 1. These for- 
mulas that derived from Sellmier equation, gives no mi- 
croscopic perspective of interactions of photon-electron 
when light passes through the material; additionally in 
sub-micron scales, RI loses its stabilization and these 
equations are not useful [15]. 

Electromagnetism have also a dispersion equation: 



for position vector of the electron, and 

P = (jPx^Py^Pz) 

for the electron momentum vector [22j. This formula 
is the simplest form in thermal equilibrium and is very 
complex in real world. This means that for attaining 
RI from quantum electrodynamics theory, we require to 
measure electron momentum vector and certainly, these 
elements can be calculated on the basis of semiclassical 
arguments. 

Now, we suggest a microscopic model for RI and wave- 
length relationship. In this pattern, average photon flux 
is related to classic electromagnetic concepts like: in- 
tensity, power and energy. For example if we have a 
monochromatic laser beam with frequency v and inten- 
sity I, then we may attain average photon flux from re- 
lation [23J: 

now if we assume thin single crystal film of NPP in 
"6" direction of crystal (or "2;" axis) radiated by a He- 
Ne laser with: A=633 nm, average power=10mw and 
beamwidth=20 microns, then from ([4]) average photon 



flux is equal to 10 photons /{s 



cm 



that signifies 
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where there are N]^ electrons per unit of volume, whose 
natural frequency is uj]^ and whose damping factor is 7/^, 

0. 



in every second 10^^ photons arrive to each centimeter 
square. Moreover from data of crystal in section (//) in 
every 36.5(A° ) in z direction, one NPP molecule exists. 
Therefore in every second 36.5 x 10^ photons interact 
with any molecule or in other words in every 27ns, one 
photon interacts with any NPP molecules. In each in- 
teraction between photon and electron in every layer of 
crystal, we suppose delay time equal to (ith layer of 
crystal). Total retardation time for m layers in crystal 
region is equal to: 



Table 1: Sellmeier data of NPP 





Ledoux et al.[3j 


selmeier form 


Datta and Banfi et al.[5j selmeier form 
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B 


C 


D 


E 


rix 


2.3532 


1.1299 


0.1678 


0.0392 


2.41704 


1.08674 


0.16933 


-0.34200 


10 


Uy 


2.8137 


0.3655 


0.2030 


-0.0816 


2.76667 


0.37156 


0.20289 


0.47880 


10 




2.1268 


0.0527 


0.1550 


-0.0608 


2.19965 


0.00457 


0.17160 


0.59363 


10 



Consequently required time for photon transmission in 
L length of crystal is equal to r, achieved from relation: 

or 

m 

n-l = |^r. (6) 

where Cq is velocity of light in vacuum. By using this 
relation, we can relate macroscopic quantity n to micro- 
scopic quantity r^. According to (3) if, n = 3,L = 3/i, 
consequently Sr is equal to 10~^^sec. Because in b 
direction of crystal in 3/i length, approximately 4024 
molecules exist, therefore the average quantity of r: 



is in order of 10~^^ sec. From this quantity and pre- 
vious quantity that in every 27ns, one photon inter- 
acts with one NPP molecule, we conclude that in every 
moment, just one photon interacts with one molecule. 
Because NPP molecule has delocalization electrons, (or 
TT-electron system), in benzene ring, that photon inter- 
acts with this electron type and it is annihilated [24J. 
We call this photon, a successful photon, (that does not 
produce phonon). By this method, we are able to sim- 
ulate microscopic perturbation of RI for NPP success- 
fully. Equations (1), (2) and (3) use average quantities 
in frequency domain, but in (6), we utilize momentary 
quantities in time domain that perturb on average val- 
ues. In our nano-quantum photonic model, we may esti- 
mate photon-electron mutual actions in very short time, 
on the order of aUosecond\ (10~^^5). We claim that 
our approach is applicable for optical material that their 
crystal structure is known, especially aromatic optical 
crystal. Because, they have a delocalization electron 
and successful photon interact with it. For inorganic 
optical material if we know active electron in optical 
phenomenon, we are able to use this method. In our ap- 
proach, we do not have any measurement and we have 
simulation and estimation with an ideal observer. This 
is the first time that we have focused on a real mate- 
rial for approximating its constitutive parameters (like 
n,£:,<j, /i). Neither in classical nor in Q.M. approach we 
only have typical material. In Q.E.D, we have not seen 
yet practical or predictable data for a real material. Our 
approach is quasi-classic and is similar to the Compton 
effect but in greater wavelengths. We have used some 
approximations for problem solution. 



4 PDF calculation of 7r-electron 
in benzene ring of NPP 

To obtain 7r-electron wavefunction for benzene molecule 
the Schrodinger equation may be solved. Since this is 
very complex, this cannot be done exactly, an approxi- 
mated procedure known as Hiickle method must be em- 
ployed. In this method, by using Hiickle Molecular- 
Orbital (HMO) calculation, a wave function is formu- 
lated that is a linear combination of the atomic orbitals 
that have overlapped [20j (see Fig. 3) (this method is 
often called the linear combination of atomic orbitals); 
that is: 

* = ^Ci$i 

i 

where the refers to atomic orbitals of carbon atoms 
in the ring and the summation is over the six C atoms. 
The Schrodinger equation for a delocalized electron is: 

Zm ^-^ 

that Vi is the atomic potential of iih atom. Following the 
common approach in quantum mechanics, we multiply 
the above equation from the left by and integrate over 
space. By solving this set of algebraic equations, the Ci 
coefficient, the wavefunction and energy E will be 
obtained, (Fig. 3). Thus, the \Ci^ is the probability of 
the TT-electron at iih atom. Thus: 

+ |C2|2 + |C3|2 + |C4|2 + |C5|2 + ICel^ = 1 
In the case of Benzene molecule: 

as followed from the symmetry of the ring [20l [25l [26] . 
But for NPP molecule there isn't such as Benzene 
molecule. NPP is a polar molecule. Nitro {NO2) is 
more powerful electronegative compound than prolinol 
and pulls 7r-electron system; consequently, the probabil- 
ity of finding TT-electron system at various carbon atoms 
of main ring isn't same and the probability of finding 
TT-electrons near the Nitro group is greater than near 
the prolinol group. Therefore there is no symmetry for 
NPP and electron cloud is spindly or oblong, (similar 
to dom-bell) (Fig. 6). We estimate this form of elec- 
tron cloud by an ellipse that our calculations would be 
uncomplicated. According to Bohr atomic model, we 
assume TT-electron is rotating around the center of posi- 
tive charge in agreement with Kepler law. The positive 




Figure 6: Assumed 7r-electron orbit of NPP molecule 
that obtained from fig. 5 for Benzene molecule (approx- 
imately). 
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Figure 8: PDF of presence probability of 7r-electron 
around the assumed orbit in fig. 7; 



Figure 7: Bohr atomic model and Kepler law for tt- 
electron system on an orbit. 

charge in this approximation is virtual. 

For attaining probability of electron presence on an 
orbit (fig. 7), we say, T time is required by radial vector 
to sweep total ii.u.v interior area of ellipse (u and v are 
semimajor and semiminor axis of ellipse respectively), in 
t times, this radial vector sweeps 

t 

TT.U.V.— 

T 

area of ellipse, (see fig. 7). If t is the time, that electron 
sweeps radian of orbit then t is obtained from this 
relation [27j : 

t = — {2arctan W- tan -f^} 

27r^ l^e ^2^^ 1 + e. cos(i9) ^ 

(7) 

Where s is ellipse eccentricity. By using this relation, we 
attain the required time for electron to traverse from to 
6-\-d6 (dt) and it is divided by total time T. By this ap- 
proach, we can determine the PDF approximately. Fig. 
8 and Fig. 9 show PDF and CDF(Cumulative Den- 
sity function) respectively for presence probability of tt- 
electron around the assumed orbit in fig. 7. Such as 
shown in this figure the PDF in apogee (near the Nitro 
group), is maximum and in perigee (near the prolinol 
group) is minimum. 

5 T Calculation 

The angle between Y vector and charge transfer action 
(N1-N2) is 58.6° and X and Z axis is perpendicular to 
Y (Fig. 10). We consider propagation along Z direction. 




Teta in Radian 

Figure 9: CDF of presence probability of 7r-electron 
around the assumed orbit in fig. 7; 




Figure 10: X, Y and Z axis and NPP molecule in dielec- 
tric frame. 




in effective center of positive charge. Consequently the 
time required for electron to go up and down is equal to: 



Vox V^Jwjn.cos{6) 
~V " KZe^ ' 



and similarly: 



that 



V2hu.m.sin{0) 



(8) 



(9) 



Figure 11: Calculation of r 



1 + e.cos{0) 

where e is the elliptical eccentricity and u is the semi- 
major axis of the ellipse. By performing Tx and Ty sepa- 
rately in every layer and using (6), we attain Ux and Uy 
separately. 



We spot a photon interacts with 7r-electron of NPP in 
first layer. After interaction, this photon gives its energy 
to electron and is annihilated. Electron absorbs energy 
and digresses in direction of photon momentum. Elec- 
tron with photon energy, may not be unbounded and 
after arriving to apogee of digression, it returns back 
to main orbit, because the photon energy is equal to 
hu = 1.96ev (h is Planck's constant and u is frequency 
of laser beam) whereas energy for ionization is greater 
than 5ev. When electron returns to main elliptical or- 
bit one photon is produced. The time coming up and 
down is r delay time. This photon after freedom goes to 
second layer in direction of annihilated photon (nonce, 
we assume the polarization doesn't change), in second 
layer this photon interacts with another delocalization 
TT-electron certainly, because the effective range of pho- 
ton is approximately equal to its wavelength and is very 
greater than the distance between molecules. This ac- 
tion is repeated for each layer. The location of photon- 
electron interaction is significant in every molecule and 
it is effective on r quantity directly. We assume that in- 
teracting photon has a field in X direction and electron 
subject to virtual positive charge center. According to 
fig. 11 we have: 



hv 



and consequently: 



1 ,.2 1 , Vox 

-mV = -mi -— 

2 2 ^cos{0) 



Vox = \l—.cos{0) 
m 



that ^mV'^ is the kinetic energy of electron after inter- 
action with photon, and columbic force, F, is equal to: 



K 



ma 



where e,a and m are charge, acceleration and mass of 
electron respectively. Ze is positive charge equivalent 



6 Attaining of 
crystal 



and rio, for NPP 



We are now in a position to take up Ux and Uy by 
Monte-Carlo method, then we simulate random num- 
ber producer using MATLAB program. This program 
produces PDF quantities explained in section (4) and 
relates each of them to every molecule. These values are 
indexing 7r-electron positions in each layer, by assump- 
tion a reference point (see fig. 7). Additionally we have 
used a MATLAB program for Monte-Carlo simulation. 
The inputs of this program are: 

1. The wavelength of incident optical beam in which 
we want to obtain refractive index of NPP crystal; 

2. /i, m, /c = Co that are Planck's constant, elec- 
tron rest mass, elementary charge. Coulomb constant 
and speed of light respectively. 

3. Unit cell parameters of NPP crystal: a, b, c, /3 and 
its other parameters that have given in section 2. 

4. L: crystal thickness that in our simulation it is 3/im. 
And the outputs of MATLAB program are: rix^ riy in 
each wavelength and their related errors from experi- 
mental data. 

System calibration is done in this method that we ob- 
tain two refractive indexes in two wavelengths with e (ec- 
centricity), u (semimajor axis of ellipse) and Z (equiv- 
alent positive charge) in a way that refractive indexes 
in two wavelengths are very close to experimental data. 
Then we would see that refractive index in other wave- 
lengths with same e, u and Z will be achieved. Of course 
these values, e, u and Z would be close to experimen- 
tal structure of crystal, for example u would be greater 
than and smaller than minimum and maximum sizes of 
six lengths of benzene hexagonal respectively, or e would 
be small and greater than zero. In other hand these val- 
ues must be logical. From this method in our simulation 
we have obtained e = 0.26, Z = 3.9, = 1.4^4^ that is 
very close to experimental and structural data . Tables 
([2]) and ([3]) compare the results of simulation with ex- 
perimental data. From these tables we perceive Ux and 



Table 2: Comparison of experimental data and simulation results for n- 



Wavelength{X — nm) 


509 


532 


546 


577 


589 


633 


644 


690 


1064 


1340 


Ux {Experimental) [3j 


2.355 


2.277 


2.231 


2.153 


2.128 


2.066 


2.055 


2.051 


1.926 


1.917 


Tlx {Simulation) 


2.290 


2.271 


2.284 


2.195 


2.200 


2.176 


2.150 


2.113 


1.910 


1.806 


Error% 


2.7 


0.23 


2.4 


1.9 


3.3 


5.3 


4.6 
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5.8 



Table 3: Comparison of experimental data and simulation results for n. 
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ny {Simulation) 
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1.913 


1.733 


1.656 


Error% 
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4.7 
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Figure 12: Comparison of experimental data, Simulation Figure 13: Comparison of experimental data, Simulation 
data and Sellmeier data for nx data and Sellmeier data for ny 



Uy that we attained from simulation, have 5.8 percent 
maximum error and 0.23 percent minimum error. These 
values are acceptable measures. Figures [12] and [13] com- 
pare these data with Sellmeier data that are provided 
from 



7 Conclusion 

The suggested physical model in this paper that is based 
on assumption of energetic particles could be a powerful 
tool for analyzing and explaining processes that happen 
in waveguides in microscopic sizes. From similar models 
we can explain optical events such as second harmonic 
generation and optical polarization in fiber optic appli- 
cations and optical communication. We will try to con- 
tinue applying this method for justification of another 
phenomenon. Our approach is quasi-quantum and its 
advantages are: 

1- We may to obtain scientific data of a real material, 

2- We utilize momentary quantities in time domain, 

3- Our method may predict some events in nanostruc- 
turs and very short time by considering short-range in- 
teratomic forces and its effect on photons. 
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